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ABSTRACT 

Context. Capella is among the few binary stars with two evolved giant components. The hotter component is a chromospherically 
active star within the Hertzsprung gap, while the cooler star is possibly helium-core burning. 

Aims. The known inclination of the orbital plane from astrometry in combination with precise radial velocities will allow very accurate 
masses to be determined for the individual Capella stars. This will constrain their evolutionary stage and possibly the role of the active 
star's magnetic field on the dynamical evolution of the binary system. 

Methods. We obtained a total of 438 high-resolution echelle spectra during the years 2007-2010 and used the measured velocities 
to recompute the orbital elements. Our double-lined orbital solution yields average residuals of 64 m S" 1 for the cool component and 
297 ms~' for the more rapidly rotating hotter component. 

Results. The semi-amplitude of the cool component is smaller by 0.045 km s _I than the orbit determination of Torres et al. from data 
taken during 1996-1999 but more precise by a factor of 5.5, while for the hotter component it is larger by 0.580 kms~' and more 
precise by a factor of 3.6. This corresponds to masses of 2.573±0.009 M Q and 2.488±0.008 M Q for the cool and hot component, 
respectively. Their relative errors of 0.34% and 0.30% are about half of the values given in Torres et al. for a combined literature- 
data solution but with absolute values different by 4% and 2% for the two components, respectively. The mass ratio of the system is 
therefore q = M A /M B = 0.9673 ± 0.0020. 

Conclusions. Our orbit is the most precise and also likely to be the most accurate ever obtained for Capella. 
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1. Introduction 

The spectrum of the Capella (a Aurigae, V = 0.07 mag) binary 
has been studied for over a century and has been observed at 
practically all wavelengths from X-rays to the radio. The liter- 
ature is too numerous to be cited here, but Torres et al. (2009 ) 
contains references to almost all of the significant work, in par- 
ticular those relevant to the present paper. The studies of Barlow 
et al. ( 119931 ) and Hummel et al. (1994 ) of all aspects of the astro- 
metric and spectroscopic orbits are cited here to emphasize their 
historical importance. 

Capella is one of the few binary systems with two moder- 
ately evolved giants, most often cited as G8 III for the primary 
and G0-1 III for the secondary. Atypical of RS CVn binaries, 
the hotter component of Capella appears to be the more active 
star rather than the cooler, usually rapidly rotating component. 
This is unsurprising for its spectral type and evolutionary status 
as this places the hot component in the Hertzsprung gap where 
rapid redistribution of the internal angular momentum is taking 
place, eventually followed by dredge up of high-angular momen- 
tum material even to the surface. Our overall goal is to learn 
more about this component by mapping its surface spots. For 
comparison, we previously investigated a very similar but single 
GO III star (31 Com; Strassmeier et al. 2010a). Quite unexpect- 
edly, we found a large polar spot in Doppler images despite the 
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predicted depth of its convective envelope of less than 20% be- 
ing much too small to account for flux-tube deflection towards 
the rotational poles. The star's rotation period was found to 
be 6.80±0.06 days from MOST data. Radial-velocity variations 
with a full amplitude of 270 ms _1 and a period of 6.76 d were 
detected from our STELLA spectra, which we also interpreted 
as being caused by the stellar rotation and made us confident 
in the use of STELLA for this high-precision work. The stellar 
mass and age for 31 Com were determined to be 2.6±0.1 M Q 
and a;540 Myr from a comparison with evolutionary tracks. To 
first order, we would expect similar values for the hot secondary 
component of Capella. 

In this paper, we report on our ongoing monitoring of 
Capella and present a new improved orbital solution. The sec- 
ondary's comparably featureless spectrum with respect to the 
cooler component and its large rotational line broadening of 
35 kms -1 make its observation difficult in spite of the favorable 
light ratio cool/hot component of «0.9 at visual wavelengths. 
Our own earlier work on Capella was based on either snap- 
shot spectroscopy (Strassmeier & Fekel 1990) or Ho- photometry 
(Strassmeier et al. 2001). The latter paper presented the detec- 
tion of a 106±3 d photometric period for the late-G giant, which 
we interpreted to be its rotation period synchronized to the or- 
bital motion. A period of 8.64±0.09 d was found for the G0-1 
giant, although it is less well defined and certainly less accurate 
than the formal error suggests but in agreement with an indepen- 
dent determination from the He i equivalent widths by Katsova & 
Scherbakov ( 1998 ). However, the knowledge of accurate abso- 
lute dimensions, in particular the component masses, radii, and 
rotational periods, are essential for comparing a system's angular 
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momentum evolution with model predictions. Its present loca- 
tion in the Hertzsprung gap suggests that there are rapid changes 
in its internal structure with a deepening convection zone and as- 
sociated changes in the total stellar moment of inertia (e.g. Kim 
& Barnes 2008). Deep mixing processes at this stage are also 
predicted to have a direct impact on the visible surface rotation 
and chemical abundance (Bohm-Vitense 2004). Precise astro- 
physical parameters of stars in this evolutionary stage may thus 
help us to understand the angular-momentum loss in late-type 
stars. 

We employed our new robotic spectroscopic facility 
STELLA in Tenerife, Canary Islands to continuously monitor 
Capella for a total of nearly four consecutive years with a ca- 
dence of one spectrum for every useable night. The instrument 
and the data are described in Sect. [2] A new SB2 orbit with, so 
far, unprecedented accuracy is presented in Sect. [3] We discuss 
its details in Sect. [4] In a forthcoming paper, we present the first 
attempt to derive a Doppler image of the surface of the rapidly 
rotating secondary component. 
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2. STELLA/SES spectroscopy in 2007-2010 

High-resolution time-series spectroscopy was obtained with 
the STELLA Echelle Spectrograph (SES) at the robotic 1.2- 
m STELLA-1 telescope in Tenerife, Spain (Strassmeier et 
al. 2010b ). A total of 438 echelle spectra were acquired over the 
course of 3 1/2 years (1312 d). The SES is a fiber-fed white-pupil 
echelle spectrograph with a fixed wavelength format of 388— 
882 nm. Despite increasing inter-order gaps in the red, it records 
the range 390-720nm continuously. Its two-pixel resolution is 
R = 55 000. The CCD is currently an E2V 42-40 2048x2048 
13.5//m-pixel device. 

Integrations on Capella were set to exposure times of 15 sec 
and 60 sec (until end of 2007) and achieved signal-to-noise ratios 
(S/N) in the range 200-500:1 per resolution element, depend- 
ing on weather conditions. The spectra were obtained between 
July 2007 (JD 2,454,305.7) and March 2011. Just 14 primary 
and 7 secondary measurements were discarded because their 
data were of too low S/N due to clouds. Our data are automati- 
cally reduced and extracted with the IRAF-based STELLA data- 
reduction pipeline (see Weber et al. 120081 1 . The two-dimensional 
(2d) images were corrected for bad pixels and cosmic rays. Bias 
levels were removed by subtracting the average overscan from 
each image followed by the subtraction of the mean of the (al- 
ready overscan-subtracted) master bias frame. The target spectra 
were flattened by dividing by the master flat, which was nor- 
malized to unity. The robot's time series also includes nightly 
and daily Th-Ar comparison-lamp exposures for wavelength cal- 
ibration and spectrograph focus monitoring. A continuous mon- 
itoring of the environmental parameters inside and outside the 
spectrograph room, most notably of temperature and barometric 
pressure, allows us to apply proper corrections. For details of the 
echelle data reduction with particular emphasis on the tempera- 
ture and pressure dependences of the SES, we refer to Weber et 
al. (120081) . 

Twenty-two radial velocity standard stars were observed 
with the same set-up and analyzed in Strassmeier et al (|2010c 
1201 H . The STELLA system appears to have a zero-point off- 
set with respect to CORAVEL (Udry et al. [19991 of 0.50 ± 
0.21 kms~', and a comparable offset of 0.46 ± 0.64 kms -1 with 
respect to the 17 stars in common with Scarfe (2010). The high- 
est external rms radial-velocity precision over the four years of 
observation was around 30 ms _1 for late-type stars with narrow 
spectral lines. All velocities in this paper are on the STELLA 



Fig. 1. Residual velocities from simulated spectra. Circles de- 
note the final setup using 62 echelle orders. Pluses indicate the 
residuals when using only 19 orders. A simulated spectrum was 
computed for the time of every observed spectrum, and the same 
method as used for the observations was applied to derive the ra- 
dial velocities of the two components. The offsets shown here are 
the differences between the velocities used to construct the simu- 
lated spectra and the values measured. Systematic offsets mostly 
affect the secondary star, with absolute values peaking at about 
50% of the single measurement error of typically 300 ms _1 
for component B. Increasing the number of orders and thus the 
wavelength range clearly suppresses the systematic errors. 

Table 1. Barycentric STELLA radial velocities of Capella. 
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zero-point scale. The individual velocities are listed in Table [T] 
available only in electronic form via CDS Strasbourg. 

3. New spectroscopic orbital elements 

3. 1 . High-precision radial velocities 

During the initial year of STELLA operations, the external 
rms values were significantly larger (120 ms ) than there- 
after (30 ms -1 ) and included a radial velocity offset of about 
300 ms~'. The final radial velocities are barycentric and cor- 
rected for Earth's rotation. 

The STELLA velocities listed here are determined from a si- 
multaneous cross-correlation of 62 echelle orders with a proper 
synthetic spectrum. To perform this, we selected from a pre- 
computed grid of synthetic spectra from ATLAS-9 atmospheres 
(Kurucz 1993 ) two spectra that match the respective target spec- 
tral classifications. For Capella A, we chose the T e fi = 5000 K 
model, and for Capella B the T e ff = 5500 K model, both with 
logg = 3.0 and solar metallicity. We combined the two tem- 
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Fig. 2. The STELLA radial velocities compared with our newly computed orbit. Filled circles denote the cool component, star A, 
open circles the hot component, star B. The two top panels show the radial velocities versus HJD and versus orbital phase, the two 
bottom panels show the residuals versus phase for the primary and the secondary, respectively. The lines are computed from our 
elements in Table [2] The horizontal dashed line in the middle panel is the systemic velocity. The residuals show a phase-dependent 
behavior in both components, but the overall rms of 64 ms _1 and 297 ms _1 for the two components fit the uncertainties of the 
individual measurements quite well. 



plate spectra to one artificial spectrum using v sin i of 4.5 km s 
and 35 kms for the A and B component, respectively, and a 
macroturbulence of 5 kms -1 . We adjusted the (wavelength de- 
pendent) intensity ratio using G8 III and GO III flux-calibrated 
spectra from Pickels ( 1998 ) multiplied with 0.955. We then ap- 
plied a series of radial velocity differences between the two 
components, and computed a two-dimensional cross-correlation 
function for each of these shifts. The highest correlation value 



in the resulting two-dimensional image corresponds to the mea- 
sured velocity of component A in one dimension, and the ve- 
locity difference of the two components in the other dimension. 
To estimate the uncertainties in this method, we evaluated a se- 
ries of 1000 Monte-Carlo variations in these two-dimensional 
images using only 63% of the original number of orders, esti- 
mating <x to be 1/1.349 of the interquartile range of the resulting 
distribution. 
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To correct for systematic radial velocity offsets induced by 
the measuring technique, we created simulated stellar spectra 
from the above templates at exactly the same wavelength spacing 
as the original observations and corrected the measured radial 
velocities with these offsets. We initially performed our analysis 
with 19 echelle orders but the final analysis was performed with 
all usable orders, which turned out to be 62 out of 82 in our case. 
Fig. [T] illustrates the systematic errors and their dependance on 
the number of orders, 19 compared to 62, and the wavelength 
range used. Doing the analysis with just 19 orders resulted in a 
TTg of about 400 ms _1 smaller than with 62 orders, even when 
taking the corrections from Fig.[T]into account. This shows that 
the systematic effects are not fully taken into account by the sim- 
ulations. 

We also correct the radial velocities of the B component for 
the difference in gravitational redshift according to Lindgren & 
Dravins (2003 ). For component A, this amounts to 132 m s _1 and 
for component B to 177 ms _1 with the masses and radii from 
Table [2] i.e. a net differential offset of 45 ms _1 . After applying 
this correction, a residual offset of the respective center-of-mass 
velocities of the two components of 59 ms _1 remained. This 
small offset could be caused by the different line-shapes of the 
two stars, e.g. produced by either convective blue- or redshifts 
(see e.g. Gray 2005 ). Therefore, we applied a total correction of 
104 ms _1 to the radial velocities of the B component. The main 
effect of this correction is a small decrease in the secondary's 
velocity amplitude. We note that all of above corrections still 
amount to changes in the final parameters of smaller than 1 <x. 

3.2. Spectroscopic orbit 

Figure [2 presents our new STELLA velocities and compares 
them with the computed velocity curves. 

We solved for the usual elements of a double-lined spectro- 
scopic binary using the general least-squares fitting algorithm 
MPFIT (Markwardt 2009 ), utilizing the reciprocal values of cr 
from Table [T] as weights. For solutions with non-zero eccen- 
tricity, we used the prescription of Danby & Burkardt (1983) 
to calculate the eccentric anomaly. The orbital period was fixed 
at 104.02173 d (Torres et al. 2009). We initially assumed a zero 
eccentricity, but the final corrected data yields a best-fit solu- 
tion at e = 0.00087 ± 0.00013 and cj = 345 ± 9.7 for both the 
SB2 solution and the solution of the primary's velocities alone. 
Using ony the secondary's velocities, we obtain a zero eccentric- 
ity. The rms of the eccentric solution is 64 ms _1 and 297 ms _1 
for the primary and secondary velocities, respectively, which is 
slightly better than the 65 ms _1 and 301 ms 4 obtained from 
the circular solution. The eccentricity is also consistent with the 
values of e = 0.0 0083 ± 0.00005 and to = 334.8 ± 4.7 that 
Torres et al. ( 120091) derived from the astrometric measurements 
of Hummel et al. (1994), and practically identical to their re- 
sult when all radial-velocity and astrometric data were combined 
(0.00087 ± 0.00021 and 324 ± 14). 



4. Discussion and conclusions 

Most prior studies of the spectroscopic orbit of Capella have 
found that the early G star is the least massive component in the 
Capella system. The only exception is an early orbit by Stmve 
& Kung (1953 ) from McDonald-Observatory spectra, as noted 
by Barlow et al. ( 119931 ). Torres et al. ((2009) collected most pre- 
viously obtained spectroscopic orbits and documented the de- 
crease in the measured semi-amplitude of the secondary with 



Table 2. Spectroscopic orbital elements of Capella. 



Parameter 


Our Value 


Torres et al. 


P (days) 




104.02173+0.00022 


rperiastron (HJD 24+) 


54 389±2.8 


47 528.514+0.016 


y (km s -1 ) 


29.9378±0.0026 


29.653+0.035 


K A (kms- 1 ) 


25.960±0.0065 


26.005+0.036 


K B (kms" 1 ) 


26.840±0.024 


26.260+0.087 


e 


0.00087±0.00013 


0.0 


co (deg) 


345±9.7 


- 


a A sin i (10 6 km) 


37.134±0.0093 




a B sin i (10 6 km) 


38.392±0.034 




M A sin 3 i (M G ) 


0.8064±0.0014 




M B sin 3 i (M Q ) 


0.7800±0.0008 




N A ,N B a 


424, 431 


504, 162* 


rms A (ms _1 ) 


64 




rms g (m s~') 


297 




mass ratio, q = M B /M A 


0.9673±0.0020 


0.9903+0.0036 


inclination, i (deg) 




137.21+0.05 


M A (M ) 


2.573±0.009 


2.466+0.018 


M B (Mo) 


2.488±0.008 


2.443+0.013 


Ra (Re) 




11.87 


Rb (Re) 




8.75 



Notes. ta) N is the number of measurements for each component. {b> A 
total number of 1015 measurements were used including astrometry. 



ever increased observational precision. The difficulty in measur- 
ing the radial velocity of the secondary (of a zero magnitude 
star!) is thereby nicely demonstrated. The earliest orbital solu- 
tions had a secondary to primary mass ratio of a;0.74, which be- 
came a;0.87 by the mid eighties in the past century, a;0.94 in the 
nineties, and 0.9903+0.0036 in the study of Torres et al. (2009 ). 
In the present paper, we verify the almost unity mass ratio but 
confine it to 0.9673+0.0020 and derive masses good to 0.3%. 

The two previously highest quality data sets - Torres et al. 
(2009) from CfA observations, and Barlow et al. (1993) from 
a combination of DAO, McDonald, and KPNO observations - 
have comparable rms for the primary of 0.46 and 0.47 kms~' 
but 0.91 and 1.38 kms -1 for the secondary, respectively. The 
CfA data set appears to be significantly more precise for the sec- 
ondary velocities than the combined Barlow et al. (1993J ve- 
locities. However, the CfA velocities are not free of systemat- 
ics. Firstly, they were obtained from single-order echelle spectra 
with a wavelength coverage of just 45 A, 117 times smaller than 
our STELLA/SES data. As shown by the authors themselves, 
this introduces phase-dependent radial velocity variations of the 
secondary star of up to +1 km s , which are produced by resid- 
ual blending and spectral lines moving in and out of the cross- 
correlation window. All CfA velocities had to be corrected for 
this effect before further use. Secondly, when the CfA velocities 
from both components were solved for separately, their best so- 
lution was found for y-velocities that differ by 0.27+0.08 km s , 
with the B-component value being the smaller one. This is not 
uncommon in double-lined solutions but this difference also per- 
sisted in their global solution. Spot activity on both stars can ac- 
count for such an effect but should statistically cancel out for a 
long enough database timeline unless phenomena such as active 
spot longitudes or variable spot lifetimes play a role. 

Our STELLA data are not affected as strongly by the above 
problems, the precision of an individual measurement even when 
compared to the best previous data being already higher by a 
factor 3^1, owing to spectrograph stability and higher spectral 
resolution. Our rms from the orbital solution is 64 ms _1 for 
the primary and 297 ms _1 for the secondary. The sampling is 



4 



M. Weber and K. G. Strassmeier: The spectroscopic orbit of Capella revisited 



also unprecedented and a total of 424 velocities for the primary 
and 431 for the secondary were available for our analysis. If we 
use the best-fit orbital period from the STELLA data alone, i.e. 
104.0243+0.0007 d, we find a practically indistinguishable orbit 
from that listed in Table|2]with rms-values of within 1 ms^'and 
4 ms^'for the cool and hot components, respectively. Because 
our time coverage is shorter than that considered in the combined 
Torres et al. orbit, 1312 d versus 5037 d, we chose to adopt their 
orbital period from their combined solution. Taking our velocity 
offset of +0.503 km s with respect to CORAVEL into account, 
our y velocity is smaller by 0.218 kms -1 than the Torres et al. 
value of 29.653+0.053 kms -1 but, more importantly, more pre- 
cise by a factor 20. The semi-amplitude of component A from 
STELLA is smaller by just 0.045 km s _1 but more precise by a 
factor of 5.5, while for component B it is larger by 0.580 km s _1 
and more precise by a factor of 3.6. As noted in Sect. 
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obtain a Kb value within the error bars similar to the one found 
by Torres et al. when using only 19 spectral orders. We there- 
fore conclude that the systematic errors in the secondary veloc- 
ity make Kb appear smaller (see Fig. [TJ. Because of the above 
rms numbers and that an inspection of the radial-velocity curves 
in Fig.[2]shows residual systematic errors only at the level of the 
accuracy of a single measurement, we conclude that our orbit is 
also more accurate. We speculate that the inclusion of the many 
different data sets in the Torres et al. solution, spectroscopic as 
well as astrometric, and their simultaneous global approach with 
a total of 26 free adjustable parameters creates high precision but 
lesser than the hoped for accuracy, possibly because of remain- 
ing systematics in the individual data sets. 

The mass ratio of a binary with two evolved giant com- 
ponents is of critical importance for the determination of their 
evolutionary status. This is particularly obvious for the Capella 
system as we have the case where two practically equally mas- 
sive giants occupy rather different locations in the H-R dia- 
gram, and yet have the same age. Mostly by virtue of its rota- 
tional line broadening, the hotter B component is the one with 
the more uncertain parameters. Its spectral classification, how- 
ever, unambiguously places it within the Hertzsprung gap where 
rapid changes in terms of surface chemistry, angular momen- 
tum etc. are taking place. The orbital elements from Table [2] 
combined with the inclination of the orbital plane from the as- 
trometric orbit of 137.21° (again cf. Torres et al.), give masses 
of 2.573+0.009 M G and 2.488+0.008 M Q for the A and B com- 
ponents, respectively. These are masses with relative errors of 
0.34% and 0.30% , about half of the errors given in Torres et al., 
and with absolute values that differ by 4% and 2% for the two 
stars. We note that we continue to denote the A component as 
the primary because it is the dominant source in the combined 
spectrum but also recall that the secondary is the brighter of the 
two components at optical wavelengths (Griffin & Griffin 1986 
Strassmeier & Fekel ll99"0i l. 
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